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AGGLOMERATION  OF  ALUMINA  PARTICLES  IN  THE  FLOW 
OF  A METALLIZED  PROPELLANT  ROCKET  NOZZLE* 


Paul  Kuentzmann 


ABSTRACT.  The  agglomeration  of  liquid  alumina 
particles  takes  place  in  the  convergent  region  and  at 
the  throat  of  the  nozzle;  it  leads  to  an  Increase  of 
performance  losses  related  to  the  velocity  and 
temperature  differences  between  the  condensed  and 
gaseous  phases. 

The  problem  is  first  schematized;  then  the 
fundamental  equations  of  the  problem  are  established 
that  give  the  evolution  of  the  distribution  law  of 
the  particle  population,  and  of  the  variables 
induced  by  the  collisions. 

Assuming  that  the  velocity  differences  may  be 
calculated  a priori,  and  choosing  a distribution 
density  in  an  analytical  form,  one  obtains  formulas 
for  the  evolution  of  the  mean  radius  and  of  the 
relative  mean  square  deviation  as  a function  of  the 
abscissa  value  along  the  nozzle  axis.  The  nozzle 
geometry  has  no  efx’ect,  and  the  agglomeration 
develops  mainly  around  the  throat. 


•Reprint  No.  7^^,  1969.  La  Recherche  Aerospatiale,  No.  131, 
July  - August,  1969.  ONERA,  Chatillon,  France.  This  study  was 
carried  out  by  the  Gunpowder  Directorate. 
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Main  Notation 


a, 

A, 


c 

A 


c 


I 
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O' 

ni 


p 

"'0 


X 


y 

p 


density  distribution  parameters; 

normal  cross  section  area  of  the  nozzle,  at 

the  throat; 

force  and  fictitious  thermal  power  induced  by 
collisions; 

characteristic  velocity; 

mass  fraction  of  the  condensed  phase; 

mass  of  one  particle; 

mans  flowrate  of  the  gaseous  phase; 

local  particle  concentration; 

chamber  pressure; 

particle  distribution  density; 

, r , radius,  average  radius,  — for  the  surface  — 
for  the  volume; 

average  temperature  of  particles,  true  temperature; 

( 0 ) 

average  velocity  of  particles,  true  velocity, 
flow  velocity  for  equilibrium; 
abscissa  taken  along  the  nozzle  axis; 
isentropic  exponent  of  the  flow  equivalent  to 
equilibrium; 

volume  mass  of  liquid  alumina; 
gas  viscosity. 


The  subscript  0 corresponds  to  initial  conditions,  the  subscript 
(?)  corresponds  to  equilibrium  flow, and  the  overlinrd  quantities 
correspond  to  flow  equivalent  to  equilibrium. 


troduction 


2y  adding  a metal  to  a solid  propellant,  the  combustion  tempera- 
ture is  increased  as  is  the  performance.  Aluminum  is  a metal  which 
is  most  widely  used  and  the  optimum  performance  is  achieved  for  a 
mass  fraction  of  the  metal  in  the  propellant  which  is  relatively 
large,  oven  though  the  flow  in  a nozzle  with  metallized  propellant 
is  the  flow  of  a biphase  mixture  in  which  the  mass  fraction  of  the 
condensed  oxide  is  on  the  order  of  0.^. 
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Capture  tests  carried  out  In  a Jet  have  shown  that  the  condensed 
phase  is  a dispersion  of  very  fine  particles  whose  radius  is  on  the 
order  of  1 y.  The  theoretical  performances  predicted  by  thermodynamic 
calculations  are  never  achieved  in  practice  and  one  can  attribute 
this  loss  partly  to  the  irreversible  exchanges  between  the  gaseous 
phases  and  the  condensed  phases.  Therefore,  in  order  to  predict 
the  performance  of  a metallized  propellant  and  to  best  design  a 
nozzle,  it  is  necessary  to  accurately  describe  the  behavior  of 
condensed  particles  in  the  flow. 

One  finds  that  the  particle  population  is  distributed  according 
to  particle  radius.  Since  it  is  easy  to  see  that  the  acceleration 
of  a particle  by  the  gaseous  flow  is  proportional  to  the  reciprocal 
of  the  square  of  its  radius,  the  particles  receive  different  veloci- 
ties and  are  subjected  to  numerous  collisions. 

Alumina  has  a temperature  close  to  3200®  K in  the  combustion 
chamber  and  is  in  the  liquid  state.  It  solidifies  around  2300°  K, 
i.e.,  in  the  divergent  part  of  the  nozzle.  The  particles  are, 
therefore,  liquid  for  quite  some  time  and  then  they  agglomerate 
due  to  their  collisions.  Their  average  mass  tends  to  increase. 

This  phenomenon  is,  therefore,  susceptible  to  reducing  the  perfor- 
mance and  it  is  important  to  determine  its  importance. 

Not  much  work  has  been  done  on  this  subject.  Marble  [l] 
established  equations  which  describe  the  variation  of  the  particle 
distribution.  A formula  is  derived  which  gives  the  growth  of  the 
average  radius  along  the  nozzle  abscissa  and  the  influence  of  the 
operational  parameters.  Crowe  and  Willoughby  [2]  numerically 
determined  the  increase  of  the  average  radius  in  the  volume  and 
studied  the  influence  of  various  parameters. 

The  equations  in  this  report  will  be  established  using  a 
different  formal  base  than  that  used  by  Marble.  The  formulas  will 
be  applied  later  on  and  we  will  establish  more  accurate  conclusions 
on  that  occasion. 
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II 


Mat henat leal  Model  — Hypotheses 


We  will  assume  that  the  particles  are  spheres  characterized  by 
their  radius.  This  hypothesis  is  not  very  realistic  because  an 
accelerated  droplet  is  deformed  because  of  the  action  of  stresses 
which  are  applied  to  its  surface  and  the  resulting  internal  circu- 
lation. Since  it  is  difficult  to  predict  the  shape  of  a droplet 
at  any  time  during  its  history,  and  therefore,  the  associated  drag, 
and  since  the  study  of  collisions  of  such  particles  would  be  parti- 
cularly difficult,  we  are  forced  to  assume  this  hypothesis.  From 
the  same  point  of  view,  we  will  neglect  the  problem  of  particle 
stability  f3]  and  we  will  assume  that  all  the  collisions  are  binary 
ones  and  only  result  in  the  production  of  a single  particle. 

We  will  assume  that  the  gaseous  flow  is  not  perturbed  locally 
by  collisions,  which  allows  one  to  calculate  the  drag  of  a particle 
independent  from  ail  of  the  others.  This  implies  that  the  average 
stay  time  of  a particle  in  the  vicinity  of  another  particle  is  very 
much  smaller  than  the  average  time  between  collisions. 

V/e  will  assume  that  the  flow  is  steady,  one-dimensional  in  any 
cross  section.  All  the  velocities  are,  therefore,  parallel  to  the 
nozzle  axis  of  revolution.  This  hypothesis  is  coherent  with  the 
collision  mechanism:  two  particles  having  velocities  parallel  to 

the  axis  will  produce  a new  particle  having  a velocity  parallel  to 
the  axis,  due  to  the  conservation  of  momentum. 

We  will  define  two  functions  of  the  abscissa  value  x and  the 
radius  r which  are  continuous;  are 

the  average  velocities  and  temperatures  of  the  particles  having  the 
same  radius.  These  quantities  are  different  from  the  true  velocities 
and  temperatures,  which  are  subjected  to  numerous  discontinuities, 
and  in  all  strictness,  depend  on  the  complete  history  of  a particle. 


These  two  latter  hypotheses  imply  that  the  Boltzmann  equation 
cannot  be  used  here.  The  integration  of  this  equation  would  be 
difficult. 

The  population  of  the  particles  will  be  locally  described  by 
the  distribution  density  f,  (r,  x),  which  is  a continuous  function 
of  r and  x.  We  assume  that  it  exists  based  on  the  large  number  of 
particles. 

Because  of  these  hypotheses,  we  are  forced  to  describe  the 
condensed  phase  as  a continuous  medium.  Our  objective  will  now  be 
to  find  the  variation  of  u^,  T^,  and  principally  f. 

III.  Variation  of  the  Distribution  Density  as  a Function  of 
Abscissa  value 

Let  us  consider  the  Op  particles  at  an  instant  t contained  in 

a unit  volume.  We  will  study  the  number  of  collisions  between  the 
particles  of  radius  a and  those  of  radius  b during  an  infinitesimal 
increase  in  time  dt. 

— The  number  of  collisions  of  one  particle  of  radius  a with 
all  the  other  particles  of  radius  b is  given  by  the  following 
elementary  formula: 


[I-, 'I  <l'' 

2 

where  n(a  + b)  is  the  collision  cross  section,  n f(b,  x)db  is  the 

P 

number  of  particles  b,  u (a)  - is  the  relative  velocity 

between  particles  a and  b. 

— The  number  of  collisions  between  the  two  particle  classes 
can  be  written  as: 
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and  all  the  collisions  produced  particles  whose  radius  is  obtained 
by  writing  down  the  equation  for  conservation  of  particle  mass: 

(a3  + 53^. 

The  number  of  particles  of  radius  r which  disappear  because  of 
the  collision  is,  therefore, 

u'\  •'.  /■  .t  ..  j ir  «■-/  , 

Jo 

II (>i  11^  Inf  ‘ iln. 

The  integral  in  the  second  term  will  be  called  I^vr,  x). 

All  the  collisions  produced  particles  of  a radius  greater  than 

r. 


Let  us  now  find  the  number  of  collisions  which  the  particles 

having  radii  between  r and  r + dr  produce  from  smaller  particles  a 

and  b.  The  integration  range  D is  found  if  we  take  the  relationship 
? ? 

a-^  + b = r into  account. 

From  this,  we  find 


... 

II  , ■ •'  '■  ' '•  . /■/>.  I II.,  . 

The  coefficient  1/2  is  introduced  in  order  to  take  into  account 
the  symmetry  property  of  a and  b in  the  integral. 

If  we  integrate  with  respect  to  B and  hold  a constant,  we  find: 

I*  «,  • »i  > ‘ 

V"'  ' |l  ''  I 


The  integral  in  the  second  term  will  be  called  x). 
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Remark 


i 

i 

db-.. 


It  is  easy  to  verify  that  the 
total  number  of  collisions  is  twice 
the-  number  of  particles  formed. 

Each  collision  produces  the  dis- 
appearance of  one  particle. 
Therefore,  we  have 


by  changing  the  variable  in  the 
second  integral. 

Let  us  now  establish  the  balance 
sheet  of  the  particles  of  the  class 
(r,  r + dr). 

At  the  time  t + dt,  we  have: 


Figure  1.  Integration  region 
for  calculating  particles 
formed  by  agglomeration. 

(dx)/[Uj^(r)]. 


tffjif.  * i!\  -r, /.  it\ 

- f/»  * i xt.  * 

particles  of  radius  r and  in  this 
expression,  we  can  replace  dt  by 


If  we  take  into  account  the  normalization  condition  for  the 
distribution  density,  we  fi.»d: 


Since 


Ur,  jr-  ,L~;  *:  / (r,  x)  -;•  ^ llr, 

we  can,  from  the  preceding  equation,  derive  (3f/3x)dx  and  tl  ^n  3f/3x 
by  letting  dx  go  to  zero.  Finally,  we  find: 


(j 


.--/I, 


(--/I,  - /*!,) 


Ur 

>~r 
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We  car.  easily  verify  the  following  relationship: 


This  equation  describes  the  evolution  of  the  distribution 
density  of  particles.  It  is  related  to  the  particle  velocity  in 
a complex  way  and  we  cannot  consider  a numerical  calculation  without 
making  it  discrete. 

IV.  Equations  for  One  Particle 

The  equation  of  motion  and  the  thermal  equation  of  a particle 
do  not  only  contain  the  classical  terms  corresponding  to  the  drag 
and  the  thermal  exchange  with  the  gas,  but  they  also  contain  an 
additional  collision  term. 

In  effect,  the  velocity  or  the  temperature  of  a particle  of 
radius  r,  which  is  the  result  of  the  agglomeration  of  two  particles 
which  are  smaller,  is  different  from  thevelocity  or  temperature  of 
the  same  particles  which  have  not  experienced  collisions.  It  is 
assumed  that  this  difference  is  instantaneously  distributed  over 
other  particles. 
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,C,;^i-'  ■.  — . .3.  ^ ...  . 


This  additional  term  must  be  zero  on  the  average,  because  we 
are  dealing  with  an  interaction. 


In  order  to  evaluate  this  additional  term,  it  is  necessary  to 
know  the  true  velocity  or  temperature  of  each  particle  before  the 
collision  and,  from  this,  to  derive  the  velocity  or  temperature  of 
the  particle  formed  after  collision.  This  true  quantity  depends  on 
the  entire  history  of  the  particle.  It  is  not  possible  to  exactly 
know  it  except  in  the  case  where  every  particle  reaches  its  asymp- 
totic regime  before  the  collision  (this  nomenclature  is  developed  in 
[3]:  hypothesis  of  ”weak  slip,"  frequently  used  by  Marble  is 

related  to  this).  If  this  is  not  the  case,  it  is  possible  to 
establish  the  existence  of  a correction  or  defect  factor  to  be  used 
for  the  true  velocity  or  the  average  velocity. 


E 


‘ 


The  momentum  of  a particle  produced  by  the  agglomeration  of  a 
particle  a and  of  a particle  b is: 


where  u’p  designates  the  true  velocity  before  collision. 

The  momentum  of  all  the  particles  of  radius  r formed  in  this 
way  per  unit  of  volume  and  during  a time  dt  is,  therefore,  given 
fay: 

-iizt'.lr.li  * | * |(<»  (/’ 

>•  Jm  'Z 

-PS. -i 

M*'"  ••  I.'''- 


The  integral  in  the  second  term  will  be  called  I^Cr,  x). 

The  variation  in  the  momentum  of  the  particles  r only  produced 
by  collisions  is,  therefore: 

• . ' Jr 
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therefore. 


5 Sf,w.;«/rrf/  (f*I,  — i*u1[r)l(r,  j:)r,). 


Since  the  resulting  force  applied  to  the  particles  r Is  given 
by: 

<ia 

lim 

4t^  tfl 

It  Is  possible  to  calculate  the  elementary  force  which  acts  on  a 
particle.  It  Is  given  by  the  following  expression: 


i*  . '*  


j-J  I|)  I. 


It  Is  easy  to  verify  that  f*p  Is  zero  on  the  average: 


that  is. 


♦*  0 ♦*  0 


which  can  be  obtained  by  changing  variables  and  from  symmetry 
considerations  In  the  first  integral. 


The  same  calculation  can  be  carried  out  In  order  to  find  the 
fictitious  elementary  power  q*p  related  to  collisions  and  relative 
to  one  particle; 


/Kir 

with 
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One  can  verify  that  q*p  Is  zero  on  the  average. 

We  will  now  use  a few  simplifications  to  find  the  variation 
of  the  distribution  density  as  a function  of  abscissa. 


These  simplifications  are  applied  to  the  velocities.  The 
principle  consists  of  assuming  that  the  particles  are  sufficiently 
small  and  that  the  flow  has  a sufficiently  small  acceleration  so 
that  we  have: 


*••1  ti  r 


..  I (hypothesis  of  weak  slip). 


One  then  shows  [3]  that  the  velocity  of  a particle  is  expressed 


by: 


•V  .J.V  7/7)- 


u^®^  is  the  fictitious  flow  velocity  for  which  kinetic  and  thermal 
equilibrium  of  the  two  phases  would  be  realized.  This  flow  corres- 
ponds to  a perfect  gas  flow  and  u^®^  does  not  depend  on  the  cross 
section  considered. 


We  will  now  write  the  velocity  differences: 

(«!  - - (fr4-=  ‘I?  «>'»  (*»- 


and  we  will  replace  Up(r)  by  u 
appears. 


(0) 


so  that  only  one  particle  velocity 


Under  these  conditions,  the  integrals  and  Ig  become: 


• V;»  rfx 

i/m  I*f  , 


J*  (r  tt)*/  (n,  x)  !<i*  — f*|  rfrt 


~ 

— (i*)4 
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As  before,  we  can  express  the  mass  flowrate  of  the  condensed 
phase  by: 


A J*  "rf  -r)  f",  (')  ",  {*) 


the  expression  for  np  becomes 


It  ^ — 

”i»  — A 


Km 


Substituting  this  In  3f/dx,  we  obtain  finally: 


2,  •< 


/(r,.r'  I,  - — (— /(^#)  I,  : 


/; 


V,  Variation  of  the  Characteristic  Parameters  of  the  Distribution 


One  distribution  law  can  be  described  by  a reduced  number  of 
parameters:  average  value,  standard  deviation,  shape  coefficients. 


By  definition,  the  average  radius  and  average  standard  deviation 
are  given  by: 

^r. 

' lo 

“ f*  I *)  dr. 


By  taking  the  derivative  under  the  sum  sign,  we  Immediately  find 


t'jr 

Hr 


at 


rfr. 


This  last  relationship  can  be  written  as: 


lijr 


i 
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a 


In  particular,  let  us  calculate  the  Integrals  In  the  second  terms: 
I ' ( - /I,  r >’1,—  /('.')  f'  l-l't  -r  ’’'li 

The  integral  In  the  numerator  Is  transformed  Into: 

/ft 

Since  we  know  that: 


we  find: 


lUMr 


I t - <1,  - f^lji  >/■  - r„-  Hr 

- - |;rVr  --f 


By  replacing  and  I^  by  their  respective  expressions,  we  obtain 

f' 

/«»  * 

II  '/•-••/  l.v-  •*  •'  I '••.  ' ' !'V  (."•  — "»  (-M  ydyh. 

.' . •• 

Dq  is  the  first  quadrant  in  the  (y,  z)  plane  and,  therefore,  this 
integral  can  also  be  written  as: 

I I lij  ;.•■ ' 7.  >1  / 'J,  »!  !»f  "•  itf  <:;! 

or  also; 


•!*  i i I.  ■"  '•*  ^ ^ ~ ** 

//»>  ' *'*  ^ ■**  ^ *’*  *'*'■  ***  ~ 

f’  j’l.  'r 
,*0 

( , *'*(«*-  “ »*)♦)•  / («•  / (('•  ~ *) 

(r*  — 
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D»  Is  between  the  r axis  and  the  first  bisector  In  the  plane 

(a,  r). 


A change  in  variables; 

y =•  (r» --«*)♦ 
: ^ a. 


tran.'if orms  this  Integral  into 


: I 


Therefore,  we  finally  have 

Ji,  ' -r^)dr 

jfj,/-’'  ' '■•/U.-r)  U',ly>- 

' Oi'*  -/-Jrfyrfc. 

Because  of  symmetry  with  respect  to  the  first  bisector,  this 
Integral  is  transformed  into: 

j J ' ■ {u„  a>  - ««.  !v } 

(ly*t  ;*)<-  (./  - .*•  - 

It  U advantageous  to  Introduce  the  following  reduced  variables: 


7.  - • 


which  finally  leads  to  these  expressions  for  the  Integral 

/<!  < » ; 


I I /i'/  <Y.  >t  / (Z.  (»,  rz.  (V,  ^i) 

;•  0^  * - I'  • y.  ~})d\<i7. 

The  calculation  is  carried  out  In  the  sane  way  for 
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Jifl- 

•'*  '■'  ^f.A*'»'* 

I’j”  (Y  i-Z)*/{Y.4./(7,.j-j(M,{Z.j»— tt,IY.jrl) 

((Y»  Z\i<  - (Y  /.I*  - 1 -■  ^)  rfVrfZ 

f'  /{Y.^)  YVY 

Let  us  now  use  the  hypothesis  of  small  slip  for  calculating  the 
velocity  differences.  We  obtain  the  following  differential  equations 

IH  * rftt*'*’  . 

•/•  rf»  " 

l-l*.  -■  Z.*  (Y* -/.=*»/ <Y.  »>/(/,  •> 

JJt.'  (,Y»  • Z»)>-(V  - Zj  . I)./Y./Z 

l‘"/(Y...  Y>.iY 


i/^*  I'"'  i i/«**'*  , 

-«•.  ...  ....  /■* 

</j  " rfr  6.V1M"'  " 

|‘j'  ,A  Z.*(Y* -Z'Sl/tV.^i/.Z.j) 


|(Y*  Z».»-  <Y'  Z».  - I ^tj-'YrfV 

f‘/(V.^)Y’«-Y 


The  calculation  of  the  Integrals,  therefore,  means  It  Is 
necessary  to  choose  an  analytic  form  for  the  distribution  density. 

VI.  First  Approximation.  Utilization  of  a Similitude  Law 

Let  us  study  laws  of  the  form 

The  normalization  condition  leads  to: 

k{t)  1 

l*r 

and  It  Is  convenient  to  select  g(X)  such  that  |„ 
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The  definition  of  the  average  radius  results  in  the  condition 

J*  Xt(X)d\  « i. 

It  is  easy  to  see  that  a is  proportional  to  r f 

in 

“ ' '"(/«  ' >)* 
and  g(X)  must  satisfy  j]^  \*^i\>jx  > t. 

Fein  [i»]  showed  that  the  particles  are  distributed  according 
to  an  exponential  law  in  the  combustion  chamber.  In  addition,  the 
results  of  captured  tests  of  particles  in  jets  carried  out  at  the 
ONERA  have  shown  the  law  to  be  of  the  form: 

Here,  we  will  assume  the  law 
defined  for  </  ■ /.  > _i. 


Prom  this  law,  we  can  derive  the  following  average  radius  and 
standard  deviation: 


^ . i-y 


such  that 


Vm  / • I*  ^ IC*/  ”* 

We  will  assume  that  b is  constant  In  this  first  approximation 
and  we  will  calculate  the  evolution  of  the  average  radius. 

Prom  this  law,  we  are  able  to  calculate  the  Integral: 
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such  that  the  differential  equation  with  respect  to  beecsaes: 

i'- ..  h'*  j t_  (*-^i)**‘‘* 

ji  •>  A;^  u*"*  i/r  r 1*  — 

f 1^,  (V  •••  Z)S(Y*  — Z*)Y*ZV'**»"V+Z» 

< (JY*  - /*)♦-  (Y  T Z)  - l)rfY«rZ. 


We  will  set 


j,  <fc> 


.'t. 


{f,  !(•*♦» 

r {b  - ifn*  i-  4) 

, (Y  ; Z}*iY*— Z*)  Y*Z*<-'*»«i‘X*Ti 


- Z»(*  — (V  Z)  i)  rfVrfZ 


In  the  preceding  differential  equation.  J^(b)  Is  a constant. 
Therefore,  the  Integration  becomes  possible.  Let  us  set: 


1 K»«  <'•'  A,  I Ju 


Jo  A !*'"*  dx 


Jt 


we  find 


Except  for  the  coefficient,  the  formula  Is  Identical  to  that 
found  by  Marble. 

Interpretation 

The  function  z Is  the  product  of  a constant  and  an  abscissa 
function. 

This  constant  can  also  be  written  as: 

1 K»<  I Kpm 
« Af  jT  B «:v 

and  characterizes  the  propellant  and  the  operational  conditions. 

The  influence  of  the  mass  fraction  of  the  condensed  phase,  the 
characteristic  velocity  for  the  gas  viscosity  (which  accelerates 
the  particles)  and  the  operational  pressure  are  clearly  demonstrated. 


The  abscissa  function  In  reality  is  a universal  function  of  the 
cross  section  (y  fixed)  and  of  the  same  u^^^  by  definition:  the 

nozzle  geometry  does  not  occur  explicitly.  In  effect,  we  have  £33: 


with 


A 

X 

(ViT-'M,  V 1 

It  is  easy  to  show  that  the  preceding  Integral  Is  limited  when 
•"  , i.e.,  when  the  equivalent  flow  Mach  number  tends  to  infinity. 

The  preceding  function  is  shown  as  a function  of  x In  Figure  2 
for  a given  nozzle  profile.  It  can  be  seen  that  the  function 
increases  rapidly  in  the  vicinity  of  the  throat  and  then  strives 
towards  a constant  value. 

Calaulation  of  J j 

By  using  polar  coordinates,  it  becomes  easier  to  calculate  the 
integral 

! I , iV  z- z--v"/V''-  »••*•»> 

(‘Y»  : z=*.»  - lY  - z.  - l).nrfZ 

I*  I » 

I I I:-' (<••14  0 • MU  0 - '•«s*ii  — ■ 

, f ♦> 

•■<►>•  fi  sin*  0 
- |{  (...,  I)  - sill  U:  - iJ.dWO 


‘-'Nozzle  profile 


Figure  2.  Universal  abscissa  function. 

The  integration  with  respect  to  R can  be  carried  out  and,  after 
simplification,  we  find: 

, 

\b  + t)  r (6  T 2)  r (*  + 4) 


f 4 (•  •t'i  0 -4-  <in  <)  fws*  0 sin*  9 
Jo  ico»  (I  — sill  Oj*  *'• 

X hib  ; f.,  iS“!«  " + 5)1  rfe. 

I cos  0 - sm  0 J 


Tranformlng  to  the  sine  of  the  double  arc,  we  find; 

, 

• " (6  -r  I)  n»  -•  211'  {b  + i) 

((5  4 + 6)  K,(6)-(6-f  5)  K,(*») 


A new  change  in  variable  finally  results  in 


J^{bl  = 


r (2  i ~ 6) 


(6r>l  l'(A  • 2)  r(/*  + -4)  2**+* 


where 


First  of  all,  we  should  note  that  can  be  calculated  for 


b = 0 


i;io)  .1 — j- 
4^4 


We  obtain  the  following  recursion  formula: 


4(6  4 i)I,'(6)~  A 

(6  i t» - — Tb.ri — ' ' 


It  is  therefore  sufficient  to  calculate  I’^  in  the  intervals 

0 and  1.  The  numerical  calculation  is  carried  out  by  expanding  the 
quantity: 


according  to  integral  powers  of  t,  which  is  possible  for  t between 
0 and  1, 


Therefore,  we  find 


- 4 ’ , ^ i (-«) 


I M —J)  {n  - Ijl  ^ 

■ n -i-  ! 


In  this  last  sum,  all  the  terras  are  negative.  We  can  verify 
that  this  series  is  convergent  using  the  D'Alembert  criterion. 
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and  It  Is  easy  to  show  that 


/lb  -\  3ilh±J)\* 

\ • l^r.:  iv  ; • 

Therefore,  It  is  necessary  to  assume  that  the  average  volume 
of  the  particles  can  decrease  for  certain  shapes  of  the  distribution 
law,  which  is  a contradiction  to  the  agglomeration  phenomenon: 
therefore,  this  analyeie  ia  inauffieient, 

VII.  Second  Approximation 


Let  us  assume  the  same  distribution  law,  but  this  time  we  will 

assume  that  b varies  with  the  abscissa  value  along  the  nozzle.  This 

2 

law  depends  now  on  two  variable  parameters  a and  b,  or  on  r^^^  and  a , 
which  are  related  to  these  two  parameters. 


The  two  equations  which  define  the*  evolution  of  the  particle 
population  are: 

i!r„  K/I/  1 

V.  iiAjt  ■ 

1 ' Z**;Y*— 

' (<Y»  rZ*)J-(Y 

. j ilr„  Kill  ^ " ti**'* 

I fir  dr  “ <>.\t  '■  u^<  IV' iTf(fr^  V) 

I fi.*  •*' 

^{Y»  --  Z*;«  — (Y*  -4  ■/*}  M -f  rf  id/.. 

Since 

o*  1 

- 6 + i - 

the  grouping  of  the  coefficient  of  b and  of  the  second  Integral 
can  be  called  Jp(b). 


We  can  introduce  the  same  var.'able  z as  before,  and  write: 


rf<j* 

'd- 


dr„ 

dz 


j,  (^'1, 


Am  ••• 

• •>*  .,  •»  J #/.l 

“ -'^1.  ,1.  - • <'»('''•• 
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By  expressing  and  o as  functions  of  a and  b,  we  find 


" =-  ■{-*  •>  ((2  i i-  :<)  »|  - {*  + 1)  Ji) 


,lh 

■(/:  ■ b, 


By  combining,  we  can  express  a as  a function  of  b: 

a 


Therefore, 


d.i  j_  (2  fc  -i  3>  J,  — (t  + 1)  J, 

lib  b + i (2 '.I  j,  — j?7  •»  ifj; 


rtb  (2  b -j  :«  J,  — (ft  + J)  J,  ^ rft  1 

" IJ  4.  (2  b',  4)  .1,  -Tfc“+  IM,  t -f  ij 


and,  consequently,  we  will  set: 

U‘t  (2  J.  --  r n J,  rf*  1 

from  which  it  follows  that: 


“‘-■"•MT.r 


By  substituting  this  expression  for  a in  the  differential 
equation  for  b,  we  find 


,11. 

rf; 


VJe  will  set 


from  which  we  find 


f*. ’-.I'.W 

■ Jo(26+ i)J,  — IlJ, 


This  relationship  defines  b as  a function  of  s and  makes  it 
possible  to  calculate  the  variations  in  the  average  radius  and  the 
mean  square  deviation  as  a function  of  abscissa  value: 


h r 1 

i.  (V 

Lli.) 

0* 

hm  *«  1 

• 

V” 

b V 1 
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Calculation  of  J 2 

The  calculation  of  3^  done  just  like  the  calculation  of  J^* 
After  transforming  to  polar  coordinates,  we  find  first  of  all 


r(26-r  6) 

* " U ~ 

n 

ft  (cm  0 — sin  6>  c<x*  6 sin*  0 

. It  fc«is  0 4-  «in  Oi**** 

((•.’  t>  - f,i  !-2  h T)  (ci!,*  0 *in*0H  — I . . . .,1 

^ i ^ .JrfO 


then 

I * 

• * (i  - 1 1 r {h  ~ i,  r“ (/n ' i»** 

f(i  b ~ «(  a' t -i  :i  / 1 1_  \ . b 

[ **“  t ; I ■ V*  i(bi  -J}  b f-lj  b'-r  ij 


and  then 


1*2  can  be  calculated  for  b « 0: 


We  can  establish  the  following  recursion  formula: 

4(*+ 

li  (*  + !»== — 36 

and  one  calculates  1*2  between  0 and  1 using  the  expansion 


rm  * ‘VI  ,,.2(2—1)...  (2-3(«-D) 


* + « + r 
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Figure  The  functions  L(b)  and  MCb). 

The  function  J2(b)  is  shown  in  Figure  3.  The  functions  L(b) 
and  M(b)  were  calculated  numerically  and  are  shown  in  Figure  4. 

numerical  Calculation 

The  variations  in  the  average  radius  and  the  average  mean 
relative  standard  deviation  were  calculated  for  the  following 
conditions: 


— biconlcal  nozzle  of  Figure  1; 

— t>Q  ■ 0 Cexponentlal  law); 


In  Figure  5*  we  can  see  that  the  average  radius  Increases 
slightly  up  to  the  geometric  throat  and  then  decreases.  The  rela- 
tive size  of  the  distribution  Increases  In  a regular  way  because 
of  the  fact  that  b decreases. 

Interpretation 

The  decrease  In  the  average  radius  can  contradict  the  Ideas 
usually  held  regarding  the  Increase  In  particles;  but  this  phenomenon 
can  be  Justified. 

a)  One  can  verify,  first  of  all ^ that  the  average  surface 
radii  and  average  volume  radii  do  Increase  and  that  this  Increase  Is 
a maximum  In  the  region  of  the  throat. 

b)  Let  us  start  with  the  expression 

i f Jo 

We  will  show  that  the  Integrals  and  I2  can  be  calculated 
when  b » 0,  l.e.,  for 


/(r,T)«  or  II  . 


In  effect,  for  I^,  we  obtain  the  following: 

I,  -■  (r  ! «•*  / (o,  r|  -u^  (rl  — u, 

r ?'««»•  I ■ (II  »)*#-»  111*— 1*1  * 

•J[i  a«-  ”,l«i 

— J*(n  ; I)»(n*~i*)e-**j 
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Figure  5.  Variation  of  the  particles  as  a function  of  the  abscissa 
value. 


and  after  calculating  the  integrals,  we  find: 

* !•  '4  djc 

IH*  - ; IP— 12  K — 2'.  T r.H»r- 9K -6)]. 


Also,  the  calculation  of  can  be  done  in  the  following  way! 


dn 


dx 


. (U»_i»,l)s<.x|.  [-  (f  -r  (R*-  Pjl)] 

fi 

’ (R*  — »»lt 


It  can  be  seen  that,  by  changing  the  variable: 

and  we  find  that  the  two  preceding  Integrals  are  equal . 


Since 


we  find 


j(/4-  (U* 


■ (K»— iV  ■' 


It 


2 « 


- f-r  . 


dx 


«*«»  — . f 

dx  • 


a , 

-»  ;tl*  • U I 'i  — f I'**  rji  j R 


By  substituting  and  I2  In  the  expression  for  3f/3x,  we 

4 

finally  arrive  at  the  relationship: 


ff  Um  I Ju 

CX  :ili  Vli  M dx 
|f  I I2K  ) 


:.i)  - n*  - ■if,  »*  ^ 2R») 

.-8i-  *»(2R*  • •«U*  * 9I»  : 6)j 


3f/3x  is  proportional  to  a certain  function  P(R)  which  Is  shown  In 
Figure  6.  We  can  see  that,  for  R « 0,  3f/3x  Is  positive.  Therefore, 
for  the  selected  analytic  form,  f(0)  cannot  have  any  other  values 
than  0(b  > 0),  l(b  * 0),  and  «(b  < 0):f(0)  Increases  starting  at  1 
and  of  necessity  we  find  that  f(0)  Is  Infinite,  l.e.,  for  negative  b. 


3f/3x  is  positive  for  the  extreme  radius  values  and  Is  negative 
for  average  values.  This  Is  to  say  that  the  proportion  of  small  and 
large  particles  tends  to  increase,  whereas  the  proportion  of  average 
particles  tends  to  decrease:  this  corresponds  very  well  to  an 

enlargement  of  the  distribution  and  allows  the  average  radius  to 
increase  or  decrease. 
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Figure  6.  Function  3f/3x  for  b « 0, 

This  aspect  of  the  phenomenon  Is  due  to  the  collision  irachanlsm 
which  depends  not  only  on  the  distribution  density  and  the  velocity 
differences,  but  also  on  the  collision  cross  section. 

c)  It  remains  to  interpret  the  results  of  capture  tests  in 
the  Jet  which  result  in  a density  distribution  corresponding  to 
b > 0. 


First  of  all,  it  is  necessary  to  know  whether  the  smallest 
particles  could  be  detected  under  a mlcr. scope  and  whether  or  not 
they  have  been  Ignored  in  the  statistics. 

One  can  also  wonder  whether  the  samples  collected  are  repre- 
sentative with  respect  to  the  real  particle  population  in  the  Jet. 

In  effect,  the  very  small  particles  follow  the  streamlines  of  the 
gaseous  flow  around  the  plate  much  more  than  the  large  ones.  The 
plate  is  located  perpendicular  to  the  nozzle  axis  and  collects  the 
particles.  The  small  particles  are  subjected  to  a greater  deflec- 
tion than  the  large  particles  and  it  would  be  appropriate  to  correct 
any  experimental  law  in  order  to  take  into  account  this  phenomenon. 
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The  simplified  calculations,  therefore,  result  In  an  adequate 
representation  of  the  agglomeration  phenomenon.  This  calculation 
attempts  to  demonstrate  the  small  role  played  by  the  nozsle  geometry 
and  to  show  the  simple  way  In  which  the  operational  parameters 
Influence  the  process.  The  conclusions  which  can  be  derived  from 
the  presented  formulas  qualitatively  agree  with  the  results  of  Crowe 
and  Willoughby  [2].  Quantitatively,  since  the  small  slip  hypothesis 
was  assumed  for  all  particles,  the  velocity  differences  and  collision 
frequencies  are  overestimated  by  the  analysis. 


VIII.  Conclusion 


The  agglomeration  of  alumina  particles  Is  therefore  primarily 

translated  by  an  expansion  of  the  distribution  law  of  the  particle 

radius  and  by  an  Increase  In  the  average  particle  surface  and  volume. 

0 

The  collisions  bring  about  a force  and  a thermal  flux  for  each 
particle  which  tends  to  reduce  the  velocity  differences  and  tempera- 
ture differences  for  particles  of  different  diameters.  The  agglomer- 
ation Is  primarily  Important  1:j  the  t.hroat  region.  Since  the  throat 
has  a great  Influence  on  the  performance  of  the  nozzle,  the  agglom- 
eration can  be  an  Important  factor  for  performance  loss.  Efforts 
must  be  made  to  utilize  the  proposed  general  formulas  In  numerical 
calculations, and  the  particle  populations  along  the  nozzle  must  be 
studied  In  experiments. 
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